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ABSTRACT 

We  present  an  approach  for  reducing  the  number  of  variables  and  constraints,  which 
is  combined  with  System  Analysis  Equations  (SAE),  for  multiobjective  optimization-based 
design.  In  order  to  develop  a  simplified  analysis  model,  the  SAE  is  computed  outside  an  op¬ 
timization  loop  and  then  approximated  for  use  by  an  optimizer.  Two  examples  are  presented 
to  demonstrate  the  approach. 
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1  Introduction 


There  are  mainly  two  classes  of  methods  for  handling  optimization-based  design  of  complex 
engineering  systems.  In  the  context  of  this  paper,  a  complex  system  entails  two  or  more 
of  the  following  features:  (i)  has  computationally  costly  SAE,  (ii)  has  a  large  number  of 
variables  and  nonlinear  constraints,  (iii)  is  multiobjective,  and  (iv)  is  decomposable  into  a 
number  of  subsystems  which  hierarchically  or  nonhierarchically  interact  with  one  another. 

The  first  class  of  methods  are  mostly  applied  to  handle  features  (i)  and  (ii)  (Vanderplaats, 
1984,  where  further  references  can  be  found),  and  more  recently  (iii)  (see,  for  example,  Zhou 
and  Tits,  1993a).  In  majority  of  the  methods  in  this  class,  the  system  analysis  is  performed 
outside  an  optimization  loop;  often  the  system  analysis  (which  is  a  part  of  an  outside  loop, 
see  Figure  2  for  an  example)  is  approximated  for  use  in  the  optimization  loop  (the  inside 
loop,  see  Figure  2)  in  an  attempt  to  reduce  the  number  of  costly  and  detailed  analyses.  The 
number  of  variables  can  also  be  reduced  by  design  variable  linking  (Vanderplaats,  1984),  this 
is  usually  done  by  problem-dependent  assumptions  (e.g.,  symmetry  in  structural  design). 
Finally,  the  number  of  constraints  can  be  reduced  by  employing  only  a  significant  subset  of 
constraints  (the  active  set )  in  the  optimization  loop  (Zhou  and  Tits,  1993b;  Gill  et  al.,  1982). 

The  second  class  of  methods  are  almost  entirely  applied  to  handle  single-objective  prob¬ 
lems  with  the  features  (i),  (ii),  and  (iv)  (  see,  for  example,  Sobieski,  1992;  Wu  and  Azarm, 
1992;  Zhao  and  Azarm,  1993,  where  further  references  can  be  found),  with  the  exception  of 
Haimes  et  al.(1990)  method  which  is  applicable  to  hierarchical  multiobjective  problems,  and 
that  of  Azarm  and  Esc.henauer  (1993)  which  does  not  handle  the  SAE  as  formulated  here.  In 
majority  of  the  methods  in  this  class,  the  system  analysis  and/or  optimization  model  can  be 
decomposed  into  a  number  of  submodels  (subsystems  or  subproblems).  The  decomposition 
(which  might  be  performed  on  the  system  analysis  or  optimization  model  or  both)  is  called 
hierarchic,  if  it  has  an  overall  tree-type  structure  with  two  or  more  levels  of  subsystems.  In 
hierarchic  decomposition,  horizontal  interaction  in  between  the  subsystems  located  at  the 
same  level  is  not  permissible.  On  the  other  hand,  the  decomposition  is  called  nonhierarchic, 
if  there  is  no  restriction  on  the  interaction  which  might  exist  in  between  the  subsystems. 

So  far,  none  of  the  methods  reported  in  the  literature  can  simultaneously  handle  all  of  the 
above-mentioned  features  which  are  becoming  increasingly  common  in  the  design  of  complex 
engineering  systems.  In  an  attempt  towards  removing  this  shortcoming,  this  paper  presents 
an  extension  of  a  recent  work  by  Azarm  and  Eschetiauer  (1993)  whereby  a  multiobjective 
approach  for  reducing  the  number  of  variables  and  constraints  is  combined  with  the  SAE 
which  might  be  hierarchically  or  nonhierarchically  coupled. 

The  remaining  sections  of  this  paper  are  organized  as  follows.  We  present,  in  Section  2, 
an  overview  of  the  definitions  and  formulation  of  the  problem.  This  is  followed  by  Section 
3,  whereby  the  solution  approach  is  presented.  The  method  is  then  demonstrated,  as  shown 
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in  Section  4,  via  two  examples:  a  simple  explosive  actuated  cylinder  (to  demonstrate  the 
solution  steps)  and  a  fairly  complex  dual-wheel  excavator.  The  paper  is  concluded  in  Section 
5  with  the  final  remarks.  Parameters  and  variables  for  the  formulation  and  examples  are 
defined  in  the  nomenclature  of  Section  8. 

2  Formulation 

As  it  was  stated  in  Azarm  and  Eschenauer  (1993),  the  overall  multiobjective  optimization 
problem  is  converted  into  a  minmax  form  in  that  the  collective  objective  is  to  minimize  the 
maximum  loss  of  all  objectives  (see  also,  Osyczka,  1984).  Hence,  the  overall  multiobjective 
optimization  formulation  is  given  as  follows: 

minmax  f(z) 

s.t.:  g(z)  <  6  (1) 

where  the  objective  and  inequality  constraint  functions,  /  and  g  components,  are  smooth 
functions  of  design  variables,  z.  The.  vector  of  design  variables,  z,  is  partitioned  into  three 
groups:  x ,  Xs ,  and  Y .  x  represents  an  yV-vec.tor  of  primary  variables,  their  impact  on  the 
design  optimization  is  assumed  to  be  global  (e.g.,  variables  that  contribute  to  many  design 
parts,  disciplines,  or  nonphysical  entities).  Xs  represents  an  5-vector  of  secondary  variables, 
their  impact  on  the  design  is  assumed  to  be  local  (e.g.,  variables  that  contribute  to  detailed 
dimensions  or  specifications  of  design  parts,  disciplines,  or  nonphysical  entities).  In  general, 
as  will  be  observed  in  the  solution  steps  (Section  3),  the  primary  variables  arc  included  into 
the  solution  process  directly,  while  the  secondary  variables  are  considered  indirectly  (i.e.,  via 
step-size  variables).  Y  represents  a  vector  of  state  variables  obtainable  as  a  solution  of  a  set 
of  simultaneous  (coupled)  equations  which  can  be  partitioned,  for  example,  into  V'n,  V  ,  Y'f, 
such  that: 

Y''  =  fa(z,Xs,Yb1Y')  (2) 

Yb  =  fb(x,X„Yc,Vn )  (3) 

Y<  =  r(i,X„Ya,Yb)  (4) 

The  above  set  of  equations  represents  the  SAE.  Each  equation,  also  referred  to  as  Contribut¬ 
ing  Analysis  (CA),  may  represent  a  particular  engineering  discipline  or  a  distinct  physical 
part  (a  subsystem)  of  a  system,  or  both  (Sobieski,  1990).  Note  that  the  division  of  SAE 
into  CAs  may:  (i)  correspond  to  physical  boundaries  present  in  the  problem  that  suggest 
its  separation  into  smaller  parts,  or  (ii)  be  purely  formal  as  in  dividing  a  set  of  equations 
into  several  subsets.  These  two  ways  of  the  SAE  partitioning  are  referred  to  as  physical  and 
nonphysical  in  the  paper. 


As  it  was  stated  before,  presumably  the  system  analysis  is  performed  outside  where  the 
objective  and  constraint  functions  (optimization  formulation)  are  evaluated.  Furthermore, 
it  is  assumed  that  in  the  optimization  formulation,  the  number  of  variables  and  constraint 
functions  (and  perhaps  the  number  of  objective  functions  as  well)  during  the  solution  pro¬ 
cess  has  to  be  reduced.  The  reduction  in  the  number  of  variables  is  achieved  through  the 
secondary  variables,  ,YS.  This  has  recently  been  developed  for  a  multiobjective  case  by 
Azarm  and  Eschenauer  (1993).  Briefly,  it  involves  partitioning  of  secondary  variables,  X, , 
into  several  groups:  Xs . ,  i  =  1  Each  XSt  represents  a  .S',-vector  of  secondary  variables 

in  an  optimization  subsystem  i  (which  might  be  an  engineering  discipline,  a  physical  part, 
or  a  nonphysical  entity).  Note  that:  S  =  ^Li  *%;  the  total  number  of  secondary  variables 
is  equal  to  the  sum  of  the  number  of  secondary  variables  in  each  optimization  subsystem: 
i  —  1,  •  •  • ,  /.  The  vector  of  secondary  variables  XSt  is  replaced  by  a  scalar  step-size  variable 
Si  as  follows: 


VHftu  Crold  i  „  j 

X3,  =  X,,  +  »i  d, 


(5) 


where  d,  (as  given  in  Section  3)  is  a  descent  direction  in  the  subspace  of  active  constraints 
which  is  obtained  from  the  optimization  subsystem  post-optimality  sensitivity  analysis.  In 
general,  the  optimization  subsystems  (or  system)  can  be  different  from  the  analysis  subsys¬ 
tems  (here  we  use  the  terms:  optimization  subsystem  and  analysis  subsystem  to  distinguish 
them  from  each  other). 

In  order  to  reduce  the  number  of  constraints,  the  cumulative  function  (or  the  KS  function) 
introduced  by  Kreisselmeier  and  Steinhauser  (1979),  which  was  also  used  by  Sobieski  et  al. 
(1985),  is  utilized  here.  As  an  example,  g^s  represents  a  cumulative  function  of  g\,-  ■  ■ , gj: 


^2exp{pgj) 


(6) 


where  p  is  a  user  controlled  coefficient.  Sobieski  et  al .( 1 988)  were  first  to  develop  a  technique, 
based  on  the  KS  function,  to  convert  a  multiobjective  optimization  problem  into  a  single 
objective  form.  The  KS  function  can  also  be  applied  to  convert  a  multiobjective  optimization 
problem  formulated  in  a  minmax  form  into  a  single  objective  form.  Figure  1  shows  an 
example  of  this  whereby  /*,  and  gks  are  the  KS  functions  for  the  objectives  and  constraints, 
respectively.  Note  that,  as  shown  in  the  Figure  1,  depending  on  the  value  of  p  there 
might  be  a  gap  in  between  the  minimum  of  /*,  and  the  minmax{f i,  /2,  /•}),  points  B  and 
A,  respectively.  These  two  points  however  will  coincide  when  p  — >  oo.  Also,  note  that  the 
choice  of  the  KS  function  as  a  means  to  reduce  the  problem  to  a  single  objective  form  implies 
that  the  objective  that  has  the  steepest  slope  after  normalization  controls  the  optimization. 


Finally,  the  SAE  can  be  approximated  in  the  optimization  formulation  in  which  the 
number  of  variables  and  constraints  are  reduced  as  well,  i.e.,  the  optimization  problem  of  (1) 
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Figure  1 :  KS  Function  for  Objectives  and  Constraints 

can  be  rewritten  as  follows  (5  and  5  are  varied  while  Y  fixed): 

minmax  f  •  •  • ,«/,  F) 

s.t.:  gks  -,3/,?)  <  0  (7) 

2:  ZCZC  Z  '  m  * 

where  Y  6  [Y^Y u],  V"  represents  a  linear  approximation  of  the  Y  which  is  computed  by 

2.  C  2.  C 

the  SAE,  Yt  and  Yu  represent  lower  and  upper  move  limits  on  the  current  value  of  Y.  The 
advantage  of  (7),  when  compared  with  (1),  is  that  the  number  of  constraints  is  reduced  to 
one  and  the  number  of  variables  is  reduced  from  (JV+5+number  of  state  variables)  in  (1)  to 
(N  +  I)  in  (7).  As  an  example,  in  a  large-scale  structural  optimization  problem,  the  reduction 
in  the  number  of  variables  and  constraints  can  (among  other  reasons)  substantially  reduce 
the  computer  storage  for  an  otherwise  very  large  jacobian  matrix  needed  by  the  optimizer. 

3  Solution  Steps 

The  solution  steps  are  summarized  as  follows  (see  also  the  flow-chart  in  Figure  2): 

•  Step  (0):  Identify  primary,  secondary,  step-size  and  state  variables.  Initialize  with 

x  —  x°,  X3l  =  =  0  and  set  k=0, 

•  Step  (1):  Solve  system  analysis  equations  for  Y',  and  perform  system  sensitivity  analysis 
(CJSE  as  described  by  Sobieski,  1990)  to  obtain  dY/di  and  dY /f)Xs. 
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Figure  2:  Flow-Chart  of  the  Solution  Steps 


Step  (2):  if  k=0  then  solve  (8)  where  x  is  varied  while  Y  is  fixed: 

minmax  f(x,  Y) 
s.t.  :  gk„(x,  ?)  <  0 

otherwise  (i.e.,  if  k  >  1),  solve  (9)  where  x  and  s  are  varied  while  Y  fixed: 

minmax  f(x,  .Sj,  •  •  • ,  .s/,  Y) 
s.t.:  gka  (i,  Y)  <  0 

”  -  2  C  Z.  C 

where  in  (8)  and  (9),  there  are  move  limits  on  Y ,  i.e.,  Y  €  (V'j,  VJ. 


Step  (3):  As  a  part  of  post  optimality  sensitivity  analysis  of  (8)  or  (9),  compute: 


and  find: 


Xliew  =  X, 

Si  S  t 


old 


+  -s,  d. 


(11) 


•  Step  (4):  Set  k=k+l,  repeat  steps  (l)-(3)  until  convergence  is  achieved. 

In  the  above-mentioned  solution  steps,  u>  and  A  are  Lagrange  multipliers  corresponding  to 
the  objective  and  constraint  functions  of  (8)  or  (9).  Furthermore,  Y  is  obtained  as  follows: 

?  =  Y+  'ZiidY/d: rt)Azi}  +  Yl{{dY/dX.i){dX„/dsi)Asi}  (12) 

i  i 

z.  c 

where  Y  is  the  current  value  of  Y.  Also: 

mt  or  *.)/«*)  =  Wf  or  «*.)/»*) local  +  £  «*/  ®r  gk.) I dY^dT / di)  (13) 

i=a,b.c 

rn  or  Sts)/a«a  =  £  {(dK7ax,)(a(/  or  w.)/ar)}(3*./a*i)  (uj 

i=n,6,c 

The  minimax  problems  of  (8)  and  (9)  are  solved  by  the  subroutine  FSQP  (Zhou  and  Tits, 
1991).  FSQP  (Feasible  Sequential  Quadratic  Programming)  is  a  set  of  Fortran  subroutines 
which  implements  algorithms  which  are  described  and  analyzed  by  Panier  and  Tits  (1993), 
Bonnans  et  al.  (1992),  and  Zhou  and  Tits  (1993a). 


4  Examples 

Two  examples  are  presented  here  to  demonstrate  the  method  developed  in  the  paper.  Both 
examples  are  selected  from  the  literature  where  they  are  formulated  as  a  single-objective 
problem.  They  are  revised  here  to  form  multiobjective  examples. 

4.1  Example  1:  Explosive  Actuated  Cylinder 

This  example  is  constructed  from  a  well-known  single-objective  optimization  problem,  a 
minimum  length  design  (/Q  of  an  explosive  actuated  cylinder  (Papalambros  and  Wilde, 
1979).  The  constraints  for  this  example  express  the  specifications  for:  kinetic  energy  (<7i), 
wall  stress  (<72),  and  geometry  constraints  (<7.3-<7s).  Its  formulation  is  revised  here  by  taking 
out  one  of  its  constraints  (piston  force)  in  the  original  formulation  and  including  it  as  an 
additional  objective  (/2)  in  the  problem.  As  shown  below,  both  objectives  have  been  scaled 
so  that  they  are  of  the  same  order  of  magnitude: 

minmax 

/l  =  (-1  +  Z'i)/ Lmax 
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h  =  [(1000^/4)^4^1]/^^ 

s.L:  (15) 

91  ■  (^w,n(l  -  7))/[1000z4^(41-^  -  i>,0"°)]  -  1  <  0 

=  Vc  +  (tt/4  )zizj 
V2  =  Vi+  (71-/4)22-35 

92  :  aeN/<ry  -  1  <  0 

af  =  (erf  -  cti a2  +  <7*)1/2 

°\  =  24(2.3  +  z\)K4  -  25) 

<r2  =  -24 


<7  3  : 

2.3/ Ana*  ~  1  <  0 

<74  : 

(21  ~f"  22)  /  Lmcix 

1  <  0 

9s  ■ 

25  <  23 

V 

0 

c*> . 

II 

-,5 

For  convenience,  the  nomenclature  for  this  example  is  given  in  the  appendix  (see  also  Pa- 
palambros  and  Wilde,  1979,  for  further  details). 


4.1.1  Solution 

Following  the  solution  steps  in  Section  3,  we  assume:  x  =  (21.24)',  X,  =  (22,2.3,25)'.  Y  = 
(Yn,YbY  =  (i>i,u2 )'.  We  then  initialize  with  x°  —  (0.4,0.34)',  X°  =  (1.2,0.82,0.67)',  .s°=0, 
and  k=0.  Next,  as  in  step  (1)  of  the  solution  steps,  we  select  and  partition  the  SAE  into 
two  nonphysical  analysis  subsystems:  subsystem  a  and  subsystem  b  (Figure  3  shows  the 
interaction  between  these  two  subsystems).  We  then  form  the  global  sensitivity  equations 


Figure  3:  Analysis  Subsystems  for  Example  1 


(Sobieski,  1990): 

(  —dv2/dv\ 

where: 


1  -dvx/dv2  \  (  dv\/dzk  \  _  (  dfa/dzk 

1  )  V  dv2/dzk  )  ~  dfb/dzk 


k  =  1,2,5 


/“  =  ur  +  (tt/4)2!2^ 


(16) 

(17) 


/  and  z 


Test  Cases 


a 


0.721 

0.724 

0.721 

0.722 

0.723 

0.724 

0 

0 

0 

0 

1.44 

1.44 

1.44 

1.44 

1 

1 

23.4 

23.4 

23.04 

23.04 

0.166 

0.166 

0.167 

0.167 

1  .45 
0.89 
22.5 
0.169 


Key: 

a:  original  double-objective  problem 
b:  original  double-objective  problem  was 

reduced  to  single-objective  by  KS  function 
c:  original  multi-constraint  problem  was 

reduced  to  single-constraint  by  KS  function 
d:  no.  of  variables  and  constraints  are  reduced 
e:  no.  of  variables  and  constraints  are  reduced 
with  linearization  on  g\ 

Table  1:  Summary  of  Results  for  Example  1 


fb  =  v,  +  (r/4)z2zl  (18) 

As  an  example,  GSE  for  z5  will  be  as  follows: 


With  k=0  at  step  (2),  we  formulate  (8)  with  p=150,  and  solve  for  updated  primary  variables: 
£=(0.75,2.3)*.  We  also  obtain  (as  a  by-product  of  FSQP):  (w1,u)2)=  (0,1)  and  A=0.2.  Next, 
we  go  to  step  (3),  whereby  (10)  we  compute  d\,  d 2,  d 3  for  X] ,  Xj,  X*,  respectively.  We  then 
express  the  secondary  variables  via  (11)  as  a  function  of  the  step-size  variable  .s  (new  values 
for  the  step-size  s  and  primary  variables  x  are  then  obtained  via  (9)  for  k  >  1).  We  now 
go  to  step  (4)  and  set  k=l  and  repeal  steps  (l)-(3).  This  iterative  process  continues  until 
there  is  no  further  improvement  in  the  value  of  objective  functions.  A  summary  of  results 
are  given  in  Table  1.  As  shown  in  the  table,  the  results  from  the  solution  steps  (test  case  e) 
compare  well  with  those  obtained  when  the  original  problem  is  solved  (test  case  a),  or  when 
the  number  of  objectives  are  reduced  to  one  (test  case  b),  or  when  the  number  of  constraints 
are  reduced  to  one  (test  case  c),  or  when  the  number  of  variables  (  including  the  step-size 
variable  .s)  are  reduced  to  three  and  the  number  of  constraints  to  one  (test  case  d). 
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4.2  Example  2:  Dual-Wheel  Excavator 

This  example  (Figures  4  and  5)  is  constructed  from  a  single-objective  optimization  problem 
(Wilson,  19921  ,'c  involves  a  minimum  weight  design  of  a  120-inch  threaded  hub-and-shaff 
assembly  f  a  dual-wheel  excavator.  There  are  9  variables  and  25  constraints  in  t  his  example. 
The  constraints  in  the  example  express  specifications  on:  hub  stress  ((/, .  </2 ) .  stress  on  shaft  at 
.arious  critical  positions  (<7:rflis)-  and  other  practical  constraints  (</ ly-y^s).  Here,  the  original 
single-objective  problem  (sum  of  the  weights  for  the  hubs  and  shaft)  has  been  converted  into 
two  objectives:  weight  of  the  shaft  (/]),  weight  of  the  hub  (/2),  with  both  objectives  scaled 
to  be  of  the  same  order  of  magnitude: 


minmax  {f\,fi} 

f\  —  (p.s'Tr/djjZ.s'pr.,  -f  Zsh=1  +  2Z&  +  2 z$zj 

+‘2zyZj  +  2Zczl  +  2z<)Zq  }  /  3300 

fi  =  (F.s’7r/4){['7(( 'i  +  «)2  —  -if)]  +  [(-7  +  B2)/2(D2u 
—  (21  +  a)2)]  4-  [(/?i  +  B2)/2[D2q  —  Df, ) }/5000 


s.t.  : 
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Figure  4:  A  Dual-Wheel  Excavator 


<t>z6  <t>z5  4>Z1  0z14>Db  4>z4  <j>z3 


Figure  5:  Shaft  of  the  Dual-Wheel  Excavator 


<7i9(2g)  ;  0.127 Ns  —  ^8  <  0 

020(2-2, 23)  :  z-2  -  0.95 Db  <  0 

<721(22,2.3)  •’  0.8 Db  —  22  <  0 

<722(2.3, 24)  :  S/Mr  +  Db/ 2  -  24/ 2  <  0 

<72.3(25, 2$)  :  25  4  0.0165  -  z6  <  0 

<724(21 ,  z2)  :  2 1  4-  0.0267  -  z2  <  0 

<725(21,25):  0.252,  -25  <0 

For  convenience,  the  nomenclature  for  this  example  is  given  in  the  appendix  (see  also  Wilson, 
1992,  for  further  details). 

4.2.1  Solution 

Following  the  solution  steps  in  Section  3,  we  assume:  x=(z\,  27)*,  .Y4l=(22,  23,  24,  28)‘, 
X,2=(25,  2(j,  29)'.  Note  that  the  secondary  variables  are  grouped  (optimization  subsystems) 
according  to  the  physical  makeup  of  the  excavator;  X3l  represents  shaft  dimensions  in  be¬ 
tween  the  two  hubs  while  X32  represents  those  outside  the  two  hubs  (Figures  4  and  5). 
State  variables  Y ,  or  the  SAE  (Figure  6),  for  this  example  are  established  so  that  they  only 
contribute  to  the  hub  stress  constraints  (<7,  and  <72).  The  hub  is  modeled  (Wilson,  1992)  as  a 


Figure  6:  Analysis  Subsystems  for  Example  2 
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ring  plate  of  linearly  varying  thickness  subject  to  concentrated  transverse  load  and  bending 
moment  at  its  outer  edge.  The  SAE  is  partitioned  into  three  nonphysical  analysis  subsys¬ 
tems:  subsystem  a  ( Yn ),  subsystem  b  (V'4),  and  subsystem  c  (Kf).  Note  that,  the  analysis 
subsystems  as  shown  in  the  Figure  6  are  different  from  the  optimization  subsystems.  The 
SAE  for  subsystem  a  represents: 


subsystem  a:  {  $  =  J  =  * 


) 


(21) 


where  n  is  the  hub’s  thickness  constant,  and  A  is  an  eigenvalue  computed  iteratively  (as  shown 
by  an  arched  arrow  on  subsystem  a  in  Figure  6)  by  the  following  characteristic  equation. 
The  characteristic,  equation  is  obtained  from  partial  differential  equations  established  for 
the  hub  based  on  the  excavation  loading  conditions  and  using  small-deflection  plate  theory 
(Conway,  1958;  Wilson,  1992): 


A4  —  2A3(n  +  2)  +  A2[(7i  -4-  2)2  +  n(l  —  /')  —  2m2] 

+A(n  +  2)  [2m2  —  n(  1  —  v)\  +  m2[m2  —  1  —  (3  +  im)(n  +  1)]  =  0 


(22) 


A  is  obtained  for  each  m  (only  11  terms  are  considered  here,  i.e.,  m=0  to  m=10)  of  the  de¬ 
flection  (w)  equation  in  subsystem  b.  As  shown  below,  subsystem  b  computes  the  deflection 
and  loads  (shear  force,  Qr ,  and  moments,  Mre)  on  the  hub: 


subsystem  b: 


f  Yb  —  w  -  E,n=o(^?Jl  +  +  O*3  +  Dr1')  cos{m0) 

Y'2b  =  Qr  =  ff(w,zx,z7,n) 

Yf  =  Mr  =  f%(w,zx,z7,n) 

Y4b  =  Mg  =  ft(w,zx,z7,n) 

Yb  =  Mre  =  ft{w,zx,z7,n) 


(23) 


where: 


(24) 


QT  =  D(r){d  /  dr]{{d2w  /  dr1)  -f  ( \/r)(dw/dr )  +  (1  /r2)(d2w/dQ2)) 

Mr  =  —D(r)[(d2w/dr2)  +  v((\  /  r)(dw  /  dr)  -f-  (l/r2)(02u>/002))] 

Me  =  —  D(r){{  \  /  r)(dw  /  dr)  +  (  \  /  r2)(d2w  /  02)  -f  v(d2w  /  dr1)) 

Mre  =  -(1  -  v) D(r)[(l /r)(d'2w / drdO)} 

and  D(r),  the  flexural  rigidity,  is  a  function  of  hub  radius,  r,  and  9  is  an  angular  coordinate 
on  the  hub. 

Finally,  subsystem  c  represents  the  boundary  conditions  which  can  be  used  to  find 
A,  B,C,  D  coefficients  for  the  above-mentioned  deflection  (re),  Wilson  (1992): 


Y{  =  A;  K/  =  B:  Y'  =  C;  Y4C  =  D 
w  =  /f(A,'fl,C,D,  A)  =  0 

subsystem  c:  ^  dw/dr  =  fi(A,  B,C,  D,  A)  =  0 

£(  forces)  =  /|(A,i?,C,D,A,n,Qr,Mr®)  =  0 
£(  moments)  =  f4(A,  B ,  C,  D ,  A,  n,  Mr)  =  0 


(25) 


Initially,  we  assumed:  x=(12.1,  8.9)',  X3,  =  (14.4,  16.6,  19.1)',  A'.j  =  (9.2,  12,  1.8)',  *,=0, 
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a 


/  and  z 


h_ 

h_ 

z\ 

z2 

z3 

£±_ 

z5 

z6 


f8 

z9 


0.80 

0.80 

12.8 

14.3 

15.2 

17.6 

7.5 

10 

5.9 

10 

1.3 


b 

0.79 

0.80 

12.7 

14.4 
15.2 

17.5 

6.9 

8.9 

5.9 
10 
1.5 


Test  Cases 


0.80 


0.81 


0.80 


0.81 


12.7 


11.9 


14.3 


13.6 


15.4 


15 


17.8 


18.6 


7.7 


9.  7 


11.8 


5.9 


10 


12 


1.4 


1.2 


0.81 

0.81 

12 

13.9 

14.8 
18.2 

9 

11.8 
6 

11.5 

1.2 


Key: 

a:  original  double-objective  problem 
b:  original  double-objective  problem  was 

reduced  to  single-objective  by  KlS  function 
c:  original  multi-constraint  problem  was 

reduced  to  single-constraint  by  KS  function 
d:  no.  of  variables  and  constraints  are  reduced 
e:  no.  of  variables  and  constraints  are  reduced 
with  linearization  on  g\  and  g2 

Table  2:  Summary  of  Results  for  Example  2 

s2=0,  and  k=0.  We  then  followed  the  solution  steps  in  Section  3  to  obtain  the  results  which 
are  summarized  in  Table  2  Finally,  it  should  be  stated  that  the  minimax  solution  reported 
here  (for  case  a  of  Table  2,  the  total  hub-and-shaft  assembly  weight  is  10,640  lb)  is  different 
from  that  reported  by  Wilson  (1992)  in  which  the  problem  is  formulated  in  a  single  objective 
form  (sum  of  the  weights  for  the  hubs  and  shaft  was  10,211  lb).  These  two  solutions  are 
essentially  two  different  Pareto  solutions  for  the  example  as  formulated  here  (Osyrzka,  1984). 


5  Concluding  Remarks 

In  this  paper  we  have  discussed  a  method  for  combining  a  reduction  technique  with  system 
analysis  equations  for  multiobjective  optimization  problems.  The  main  characteristics  of  the 
method  presented  are  that:  (i)  the  number  of  variables  and  constraints  can  be  reduced,  (ii) 
th*3  SAE,  which  might  be  costly  and  hierarchically  or  nonhierarchically  decomposed,  is  per¬ 
formed  outside  the  optimization  loop,  and  (iii)  it  is  applicable  to  multiobjective  optimization 
problems. 

The  method  has  been  demonstrated  by  two  examples:  (i)  a  simple  explosive  actuated 


cylinder,  and  (ii)  a  fairly  complex  dual- wheel  excavator.  It  has  been  shown  that,  for  both 
examples,  when  the  reduction  measures  (test  cases  (b)-(e)  for  both  examples)  are  applied 
they  can  usually  obtain  a  solution  fairly  close  to  that  of  the  original  problem  (test  case 
(a)).  The  small  difference  in  the  solution  could  be  eliminated,  for  example,  by  increasing  the 
number  of  variables,  constraints,  etc.,  to  resolve  a  lack  of  sufficient  degrees- of- freedom . 
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Formulation: 


di 

I 

fks 


9 

9ks 

I 

k 

N 


S 

Si 

x 

X . 

V 

ya 

2  C 

Y 


\ 


—  Contributing  analysis  equations  in  analysis  subsystem  a 
=  Descent  direction  in  optimization  subsystem  i 
=  Vector  of  objective  functions1 
=  Cumulative  objective  function 
=  Vector  of  inequality  constraint  functions 
=  Cumulative  constraint  function 
=  Number  of  Optimization  subsystems 
=  Iteration  counter 
=  Number  of  elements  in  tildex 
=  Step-size  variable  in  optimization  subsystem  i 
=  Number  of  elements  in  X, 

=  Number  of  elements  in  Xa< 

=  Vector  of  primary  variables 
=  Vector  of  secondary  variables 

=  Vector  of  secondary  variables  in  optimization  subsystem  ?’ 
=  Vector  of  state  variables 

=  Vector  of  state  variables  in  analysis  subsystem  a 
=  Lower  and  upper  move  limits  on  the  current  value  of  V 
=  Linear  approximation  to  Y 
=  Vector  of  overall  design  variables 
=  Lagrange  multipliers  for  constraint  functions 
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p 


LO 

Example  1: 


D„ 


N 


V\ ,  V-i 
Ve 

W 

nun 

~1 

Zi 

~3 


Z\ 

-5 

7 

Cy 

°\ ,  <t2 


=  User  defined  coefficient  for  KS  function 
=  Lagrange  multipliers  for  objective  functions 

=  Maximum  allowable  cylinder  outside  diameter,  1  in 
=  Maximum  piston  force,  700  lb 
=  Maximum  cylinder  total  length,  2  in 
=  Safety  factor,  3 

=  Initial,  final  volume  of  combustion,  in3 
=  Fixed  chamber  volume,  0.084  in3 
=  Minimum  kinetic  energy,  600  lb-in 
=  Unswept  cylinder  length,  in 
=  Working  stroke  of  piston,  in 
=  Outside  diameter  of  cylinder,  in 
=  Initial  pressure  of  combustion,  ksi 
=  Piston  diameter,  in 
=  Ratio  of  specific  heats,  1.2 
=  Equivalent  stress,  ksi 
=  Yield  strength,  125  ksi 
=  Principal  stress,  ksi 


Example  2: 

a  =  Diameter  constant  for  the  hub,  in 

Bx ,Bi  =  Hub  dimensions,  in 

Dio  =  Hub  dimension,  in 

Dn  =  Hub  dimension,  in 

Db  =  See  Figure  5,  in 

f,  =  Safety  factor 

Ng  —  Number  of  seals 

S/Mr  =  Radial  height  of  the  shoulder  at  the  bearings,  in 

2,  =  Shaft  dimension  (i  =  1, •  •  •  ,9),  see  Figure  5,  in 

Zb  =  See  Figure  5,  in 

Zc  =  See  Figure  5,  in 

Zsb  =  See  Figure  5,  in 

Zsp  =  See  Figure  5,  in 

A  =  Eigenvalue 

v  =  Poisson’s  ratio 

ps  —  Density,  lbm/in3 

aab  =  Alternating  stress  on  shaft  at  bearings,  ksi 

Va.end  =  Alternating  compressive  stress  on  shaft  at  end  connections,  ksi 

a atnt  =  Alternating  tensile  stress  on  shaft  at  end  connections,  ksi 

a nn  =  Alternating  shear  stress  on  shaft  at  end  connections,  ksi 

aah  =  Alternating  stress  on  hub,  ksi 

aaTev  =  Alternating  stress  on  shaft  due  to  reverse  loading,  ksi 

Vathd  =  Alternating  stress  on  ACME  threads,  ksi 

Can  =  Alternating  stress  on  shaft  at  ACME  thread  undercut,  ksi 
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&  ax  th 

ae 

&?nb 

&mend 

&ment 

&mes 

&mk 

&mrev 

&mthd 

®mu 

&mxth 

&ut 


=  Alternating  shear  stress  on  ACME  threads,  ksi 
=  Endurance  limit,  ksi 
=  Mean  stress  on  shaft  at  bearings,  ksi 
=  Mean  compressive  stress  on  shaft  at  end  connections,  ksi 
=  Mean  tensile  stress  on  shaft  at  end  connections,  ksi 
=  Mean  shear  stress  on  shaft  at  end  connections,  ksi 
=  Mean  stress  on  hub,  ksi 

=  Mean  stress  on  shaft  due  to  reverse  loading,  ksi 
=  Mean  stress  on  ACME  threads,  ksi 
=  Mean  stress  on  shaft  at  ACME  thread  undercut,  ksi 
=  Mean  shear  stress  on  ACME  threads,  ksi 
=  Ultimate  tensile  strength,  ksi 
=  Yield  strength,  ksi 
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